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Abstract 



We describe the creation and evolution of electron-positron pairs in a strong electric 
field as well as the pairs annihilation into photons. The formalism is based on 
Q"^ ' generalized Vlasov equations, which are numerically integrated. We recover previous 

■ results about the oscillations of the charges, discuss the electric field screening and 
\ the relaxation of the system to a thermal equilibrium configuration. The timescale 

O ' of the thermalization is estimated to be ~ 10^ — 10^ h/rUeC^- 

p 
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^ ' Three different earth-bound experiments and one astrophysical observation 

have been proposed for identifying the polarization of the electronic vacuum 
^ I due to a supercritical electric field [S > Sc = mlc^/eh, where rrie and e are the 

■ electron mass and charge) postulated by Sauter-Heisenberg-Euler-Schwinger 

[1]: 

(1) In central collisions of heavy ions near the Coulomb barrier, as first pro- 
posed in [2,3] (see also [4,5,6]). Despite some apparently encouraging 
results [7], such efforts have failed so far due to the small contact time of 
the colliding ions [8,9,10,11,12]. Typically the electromagnetic energy in- 
volved in the collisions of heavy ions with impact parameter li ~ 10~^^cm 
is El ~ 10~^erg and the lifetime of the diatomic system is ti ~ 10~^^s. 

(2) In collisions of an electron beam with optical laser pulses: a signal of 
positrons above background has been observed in collisions of a 46.6 GeV 
electron beam with terawatt pulses of optical laser in an experiment at 
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the Final Focus Test Beam at SLAC [13]; it is not clear if this experimen- 
tal result is an evidence for the vacuum polarization phenomenon. The 
energy of the laser pulses was E2 ~ lO^erg, concentrated in a space-time 
region of spacial linear extension (focal length) I2 ~ 10~^cm and temporal 
extension (pulse duration) t2 ~ 10~^^s [13]. 
(3) At the focus of an X-ray free electron laser (XFEL) (sec [14,15,16] and 
references therein). Proposals for this experiment exist at the TESLA 
colhder at DESY and at the LCLS facility at SLAC [14]. Typically the 
electromagnetic energy at the focus of an XFEL can be E3 ~ lO^erg, con- 
centrated in a space-time region of spacial linear extension (spot radius) 
I3 ~ 10~^cm and temporal extension (coherent spike length) ~ 10~^^s 
[14]. 

and from astrophysics 

(1) around an electromagnetic black hole (EMBH) [17,18,19], giving rise to 
the observed phenomenon of gamma-ray bursts (GRB) [20,21,22,23]. The 
electromagnetic energy of an EMBH of mass M ~ IOMq and charge 
Q ~ O.IM/ \/G is E4 ~ lO^^erg and it is deposited in a space-time region 
of spacial hnear extension I4 ~ lO^cm [18,24] and temporal extension 
(collapse time) ^4 ~ 10~^s [25]. 

In addition to their marked quantitative difference in testing the same basic 
physical phenomenon, there is a very important conceptual difference among 
these processes: the first three occur in a transparency condition in which 
the created electron-positron pairs and, possibly, photons freely propagate to 
infinity, while the one in the EMBH occurs in an opacity condition [26]. Under 
the opacity condition a thermalization effect occurs and a final equipartition 
between the e^e^ and 7 is reached. Far from being just an academic issue, 
this process and its characteristic timescale is of the greatest importance in 
physics and astrophysics. It has been shown by a numerical simulation done 
in Livermore and an analytic work done in Rome [26], that, as soon as the 
thermalization of e+e^ and 7 created around an EMBH has been reached, 
the plasma self propells outwards and this process is at the very heart of the 
gamma-ray burst (GRB) phenomenon. A critical step was missing up to now: 
how to bridge the gap between the creation of pairs in the supercritical field 
of the EMBH and the thermalization of the system to a plasma configuration. 
This letter reports some progress on this topic with special attention to the 
timescale needed for the thermalization of the newly created e^e~ pairs in the 
background field. The comparison of the thermalization timescale to the one 
of gravitational collapse, which occurs on general relativistic timescale, is at 
the very ground of the comprehension of GRBs [25]. 

The evolution of a system of particle-antiparticle pairs created by the Schwinger 
process has been often described by a transport Vlasov equation (see, for ex- 
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ample, [27,28]). More recently it has been showed that such an equation can 
be derived from quantum field theory [29,30,31]. In the homogeneous case, the 
equations have been numerically integrated taking into account the back re- 
action on the external electric field [32,33,34,35]. In many papers (see [36] and 
references therein) a phenomenological term describing equilibrating collisions 
is introduced in the transport equation which is parameterized by an effective 
relaxation time r. In [36] one further step is taken by allowing time variabil- 
ity of r; the ignorance on the collision term is then parameterized by a free 
dimensionless constant. The introduction of a relaxation time corresponds to 
the assumption that the system rapidly evolves towards thermal equilibrium. 
In this paper we focus on the evolution of a system of e+e" pairs, explic- 
itly taking into account the scattering processes e+e" ^ 77. Since we are 
mainly interested in a system in which the electric field varies on macroscopic 
length scale {I ~ lO^cm, above), we can fimit ourselves to a homogeneous 
electric field. Also, we will use transport equations for electrons, positrons 
and photons, with collision terms, coupled to Maxwell equations. There is no 
free parameter here: the collision terms can be exactly computed, since the 
QED cross sections are known. Starting from a regime which is far from ther- 
mal equilibrium, we find that collisions do not prevent plasma oscillations in 
the initial phase of the evolution and analyse the issue of the timescale of 
the approach to a e~^e~j plasma equilibrium configuration, which is the most 
relevant quantity in the process of gravitational collapse [25]. 



The motion of positrons (electrons) is the resultant of three contributions: the 
pair creation, the electric acceleration and the annihilation damping. The ho- 
mogeneous system consisting of electric field, electrons, positrons and photons 
can be described by the equations 



9t/e + eE9p/e = 5(E,p)- 



1 



(2^) 



^ep^Ce(t,p) , 



dtE = -jp 



E)-jc (t), 



(1) 

(2) 
(3) 



where /g = /e {t, p) is the distribution function in the phase-space of positrons 
(electrons), = /-y (t, k) is the distribution function in the phase-space of 

1 /2 

photons, E is the electric field, Cp = (p ■ p + ml) is the energy of an electron 
of 3-momentum p (mg is the mass of the electron) and Ck = (k • k) 



1/2 



IS 



the energy of a photon of 3-momentum k. /g and L are normalized so that 



/ fe (t, p) = He (t), J (t, k) = (t) , whcre He and are number 



densities of positrons (electrons) and photons, respectively. The term 



5(E,p) = (27r)^,^ = -|eE|log 



1 — exp 



7r(m^+pj_) 
|eE| 



(4) 



is the Schwinger source for pair creation (see [32,33]): py and pj^ are the 
components of the 3-momentum p parallel and orthogonal to E. We assume 
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that the pairs are produced at rest in the direction parallel to the electric field 
[32,33]. We also have, in Eqs. (1), (2) and (3), 

J PI •'l ^2 

X \Mf [/, (p) /e (Pl) - f, (ki) /, (k^)] , (5) 

C, {t, k) ~ J {p,+p,-k- h) 

X \M\' [fe (Pl) fe (P2) - A (k) f, (ki)] , (6) 

which describe probability rates for pair creation by photons and pair anni- 
hilation into photons, M. — A4e+(pi)e-(p2)?i7(k)7(ki) being the matrix element 
for the process e+ (pi) e" (P2) — > 7 (k) 7 (ki). Note that the coUisional terms 
(5) and (6) are either unapplicable or negligible in the case of the above three 
earth-bound experiments where the created pairs do not originate a dense 
plasma. They have been correctly neglected in previous works (see e. g. [16]). 
CoUisional terms have also been considered in the different physical context of 
vacuum polarization by strong chromoelectric fields. Unlike the present QED 
case, where expressions for the cross sections are known exactly, in the QCD 
case the cross sections are yet imknown and such coUisional terms are of a 
phenomcnological type and usefull uniquely near the equilibrium regime [36]. 
Finally j, (E) = 2f / ^e^S (E, p) and je (t) = 2en, J ^fje (p) are po- 
larization and conduction current respectively (see [28]). In Eqs. (5) and (6) 
we neglect, as a first approximation, Pauli blocking and Bose enhancement 
(see e.g. [33]). By suitably integrating (1) and (2) over the phase spaces of 
positrons (electrons) and photons, we find the following exact equations for 
mean values: 

|n, = 5(E)-rz^(ai^')e + <(^2^'%, 
l^e (ep>e = erieE • (v)^ + |E • - {e^'^iv"), + < {e^T2v")^ , 

|ne (p)^ = erieE - nl {y>(Tiv')^ , 

|E = -2ene(v),-j,(E), (7) 

where, for any function of the momenta 

{F (Pl, p„)), = nr j^-^F (Pl, pO • fe (Pl) • ... • /e (Pn) , (8) 

{G (ki, k,)), ^ < / ^-^^ G (ki, k,) • /, (kO •...•/, (k,) . (9) 

Furthermore v' is the relative velocity between electrons and positrons, v" is 
the relative velocity between photons, ai — ai (e^"^^ is the total cross section 

for the process e+e~ — > 77 and (J2 = (72 is the total cross section for 
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the process 77 e+e (here e^oM -g energy of a particle in the reference 
frame of the center of mass) . 

In order to evaluate the mean values in system (7) we need some further 
hypotheses on the distribution functions. Let us define p\\, Cp and such 
that 



P\\ 



(ep)g = ep = (p| + pi + miy/^. We assume 

/e {t, p) OC He {t) S - P\\) S (pi - pi) . 



(10) 



Since in the scattering e'^e 77 the coincidence of the scattering direction 
with the incidence direction is statistically favored, we also assume 



(t, k) OC (t) 5 (ki - ki) [6 (kii - kii) + 5 (kii + ki\) 



fill 



where k\\ and k_L have analogous meaning as p\\ and p^ and the terms 5 (^k\\ — k\\j 

and 6 (jv\\ + fcy) account for the probability of producing, respectively, for- 
wardly scattered and backwardly scattered photons. Since the Schwinger source 
term (4) implies that the positrons (electrons) have initially fixed py, p\\ — 0, 
assumption (10) ((11)) means that the distribution of (fcy) does not spread 
too much with time and, analogously, that the distribution of energies is suf- 
ficiently peaked to be describable by a 5— function. The dependence on the 
momentum of the distribution functions has been discussed in [33,30]. Ap- 
proximations (10), (11) reduce Eqs. (7) to a system of ordinary differential 
equations. In average, since the inertial reference frame we fix coincides with 
the center of mass frame for the processes e+e^ ^ 
species. Substituting (10) and (11) into (7) we find 



77, e 



CoM 



~ e for each 



S (E) - 2nlaip-^ TTey +2riia2, 



-1 



- 4:nla2, 



2n^p^a2, 



An^Pya2, 



f^TTeW = eUeS - 2ne7le\\(JlPe ^ 



dt 



£ = -2enePe ^ ^e\\ - jp {£) , 



(12) 



where Pe 



n-j,ek, TTeii = neP|| are the energy density of positrons 



(electrons), the energy density of photons and the density of "parallel mo- 
mentum" of positrons (electrons), £ is the electric field strength and jp the 
unique component of jp parallel to E. ai and (T2 are evaluated at e^oM _ ^ 
for each species. Note that Eqs. (12) are "classical" in the sense that the only 
quantum information is encoded in the terms describing pair creation and scat- 
tering probabilities. Eqs. (12) are consistent with energy density conservation: 

I (Pe + P,+ '2£') = 0. 



5 



The initial conditions for Eqs. (12) are Ue = n^, = pe = P7 = T^e\\ = 0, S = £q. 
In Fig. 1 the resuhs of the numerical integration for £q = 9Sc is showed. 
The integration stops at t = 150 tc (where tq = h/rrieC^). Each variable 
is represented in units of rrie and Ac = h/rrieC. The numerical integration 
confirms [32,33] that the system undergoes plasma oscillations: a) the electric 
field oscillates with decreasing amplitude rather than abruptly reaching the 
equilibrium value; b) electrons and positrons oscillates in the electric field 
direction, reaching ultrarelativistic velocities; c) the role of the e"'"e~ 77 
scatterings is marginal in the early time of the evolution, the electrons are too 
extremely relativistic and consequently the density of photons builds up very 
slowly (see. details in Fig. 1). 

At late times the system is expected to relax to a plasma configuration of 
thermal equilibrium and assumptions (10) and (11) have to be generalized to 
take into account quantum spreading of the distribution functions. It is never- 
theless interesting to look at the solutions of Eqs. (12) in this regime. In Fig. 2 
we plot the numerical solution of Eqs. (12) but the integration extends here 
all the way up to t = 7000 tc (the time scale of oscillations is not resolved 
in these plots). It is interesting that the leading term recovers the expected 
asymptotic behaviour: a) the electric field is screened to about the critical 
value: £ ~ £c for t ~ 10^ — IO^tq ^ tc; b) the initial electromagnetic energy 
density is distributed over electron-positron pairs and photons, indicating en- 
ergy equipartition; c) photons and electron-positron pairs number densities 
are asymptotically comparable, indicating number equipartition. At such late 
times a regime of thermalized electrons-positrons-photons plasma begins and 
the system is describable by hydrodynamic equations [25,26]. 

We provided a very simple formalism apt to describe simultaneously the cre- 
ation of electron-positron pairs by a strong electric field S > Sc and the pairs 
annihilation into photons. As discussed in literature, we find plasma oscilla- 
tions. In particular the collisions do not prevent such a feature. This is because 
the momentum of electrons (positrons) is very high, therefore the cross section 
for the process e^e~ — » 77 is small and the annihilation into photons is neg- 
ligible in the very first phase of the evolution. As a result, the system takes 
some time {t ~ 10^ — IO^tc) to thermalize to a e"'"e~7 plasma equilibrium 
configuration. We finally remark that, at least in the case of electromagnetic 
Schwinger mechanism, the picture could be quite different from the one pre- 
viously depicted in literature, where the system is assumed to thermalize in a 
very short time (see [36] and references therein). 

It is conceivable that in the race to first identify the vacuum polarization pro- 
cess a la Sauter-Euler-Heisenberg-Schwinger, the astrophysical observations 
will reach a positive result before earth-bound experiments, much like in the 
case of the discovery of lines in the Sun chromosphere by J. N. Lockyer in 
1869, later identified with the Hehum spectral fines by W. Ramsay in 1895 
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[37]. 



We are grateful to A. Ringwald and P. Chen for clarifications on the earth- 
bound experiments and general discussions, an anonimous referee for quoting 
the corresponding analysis in QCD and P. Giannone for the quotation relative 
to the Helium lines. 
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Fig. 1. Plasma oscillations. We set £q = 9£ci t < 150rc and plot: a) electromagnetic 
field strength; b) electrons energy density; c) electrons number density; d) photons 
energy density; e) photons number density as functions of time. 
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Fig. 2. Plasma oscillations. We set Eq = 9fc) t < TOOOrc and plot: a) electromagnetic 
field strength; b) electrons energy density; c) electrons number density; d) photons 
energy density; e) photons number density as functions of time - the oscillation 
period is not resolved in these plots. The model used should have a breakdown 
at a time much earlier than 7000rc and therefore this plot contains no more than 
qualitative informations. 
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